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OPERATING  RANGES  OF  METEOROLOGICAL  WIND  TUNNELS 
FOR  THE  SIMULATION  OF  CONVECTIVE  BOUNDARY  LAYER  PHENOMENA 


R.N.  Meroney^  and  W.H.  Melbourne^ 

Department  of  Civit  Enoinegf ino,  Colorado  State  University 
Fort  CoUins.  Colorado.  80S23.  U.S.A. 

^Department  of  Hechanieat  gngineering.  Monash  University 
Ctavton.  Victoria.  3168.  Australia 

Abstract:  The  operating  ranges  of  meteorological  wind  tunnels  for  convective  boundary  layer  (CBL)  simulation 
are  defined  in  this  paper  based  on  a  review  of  the  theoretical  and  practical  limitations  of  the  flow  phenomena 
and  the  facilities  available.  Wind-tunnel  operating  ranges  are  limited  by  the  dimensions  of  the  simulated 
circulations  and  of  the  tunnel  itself,  the  tunnel  flow  speed  and  turbulence  processes,  and  the  characteristics 
of  the  measurement  instrunentation.  When  it  is  desired  to  simulate  both  the  convective  boundary  layer  and  the 
behavior  of  other  flows  imbedded  within  the  boundary  layer,  such  as  power-plant  plume  rise  and  dispersion,  then 
additional  constraints  exist  on  the  fluid  modeling  process.  The  capabilities  of  meteorological  wind  tunnels 
can  also  be  extended  through  the  judicious  use  of  boundary  and  side  wall  flow  controls. 


1.  Introduction 

Convective  boundary  layer  (CBL)  circulations  are  surface  temperature- forced 
atmospheric  circulations  observed  almost  daily  over  most  of  the  earth's  surface. 
An  understanding  of  the  convective  boundary  layer  (CBL)  has  only  begun  to  emerge 
from  laboratory,  field  and  numerical  research  in  the  last  twenty  years.  Recent 
studies  have  focused  on  the  characteristics  of  the  large  and  central  mixed  layer, 
but  still  unresolved  are  the  details  of  how  the  various  CBL  sublayers  are  linked 
together.  Much  of  the  remaining  difficulties  relate  to  the  fact  that  definitive 
measurements  in  regions  above  the  meteorological  tower  heights  but  below  the 
mixed  layer  and  above  the  mixed  layer  within  the  capping  entrainment  layer  are 
difficult  to  acquire. 

Much  of  today's  understanding  of  the  CBL  has  come  from  convection  tank  and 
laboratory  studies.  Conventional  type  and  size  water  tank  and  wind  tunnel 
facilities  provided  data  which  were  used  to  interpret  the  mixed-layer  atmospheric 
physics  and  to  validate  analytic  and  numerical  models  for  plume  behavior. 
Routine  engineering  investigations  of  CBL  impact  on  power-plant  plumes  or  other 
activities  of  man  will  be  constrained  by  the  inherent  limitations  of  the  size  and 
type  of  simulations  used.  This  paper  examines  some  of  the  limitations  of  fluid 
modeling  and  critiques  the  potential  for  boundary- 1 ayer  control  methods  to  extend 
the  capabilities  of  existing  meteorological  wind  tunnels. 

a.  Global  Characteristics  of  the  CBL 

Predictions  of  many  features  of  neutral  and  stable  atmospheric  boundary  layers 
were  made  possible  through  Lagrangian  similarity  assumptions  (Monin  and  Obukhov, 
1954),  K-modeling  and  mixing  length  numerical  approaches,  and  even  laboratory 
simulation  in  air  and  water  facilities  (Arya  and  Plate,  1969).  Initially  it  was 
hoped  that  the  Monin-Obukhov  similarity  approach  would  also  explain  the  vagaries 
of  the  unstable  boundary  layer,  but,  as  height  increased  above  ground  level,  data 
departed  significantly  from  predictions.  Then  Willis  and  Deardorff  (1974) 
performed  a  series  of  experiments  on  the  behavior  of  stably  stratified  water 
layers  heated  from  below  in  a  water  tank,  and  they  found  that  the  characteristic 
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length,  velocity  and  temperature  scales  were  the  mixing  layer  depth,  z. ,  the 
convective  velocity,  w, ,  and  the  convective  temperature,  9.,  respectively;  wnere 


w* 


'  gQz,- 
pZ  T 


1/3 


and  9,  = 


Earlier,  based  on  sub-grid  scale  modeling  of  unstable  channel  flow,  Deardorff 
(1970,  1972)  had  suggested  this  possibility  existed  for  situations  where  -z^L^ 
>  10. 


Caughey  (1984)  and  Kader  and  Yaglom  (1990)  suggest  that  the  convective  boundary 
layer  may  be  idealized  as  a  multilayer  structure.  Proceeding  upwards  from  the 
surface  one  can  identify  1)  the  dynamic  sublayer  in  which  shear  plays  a  dominant 
role  such  that  all  velocities  and  lengths  scale  with  friction  velocity,  u, ,  and 
position,  z.  This  layer  persists  as  long  as  the  production  of  shear  produced 
vertical  fluctuation  energy  is  less  than  buoyant  production  of  such  motions  [z 
<  L.  =  (xu, )^/(Q4)  ] .  Next  2)  the  dvnamic-convective  sublayer  denotes  a  region 
where  horizontal  velocities  are  scaled  by  the  characteristic  shear  velocity,  u„ 
but  vertical  motions  scale  by  the  convective  velocity,  w, .  The  upper  edge  of  the 
dynamic-convective  sublayer  is  at  a  height  z  of  the  order  L  =*  /cL^ .  When  z  >  L., 
=  (k,  u,  )^/(Q)3)  then  shear  stress,  ,  is  no  longer  important  but  height  continues 
to  be  the  significant  length  scale.  In  this  3)  free  convection  sublayer  the 
characteristic  velocity  for  all  directions  becomes  the  convective  velocity,  w,. 
Somewhere  around  z  =  O.lz,-  the  free  convection  sublayer  blends  into  the  4)  mixed 
layer  for  which  the  relevant  length,  velocity,  and  temperature  scales  are  now 
mixed  layer  depth,  z,  ,  convective  velocity,  w,,  and  convective  temperature,  9,. 
Finally  between  0.8z^  and  1.2z^  exists  the  5)  entrainment  interfacial  layer 
which  depends  upon  the  convective  velocity,  w, ,  the  depth  of  the  capping 
inversion,  and  the  strength  of  the  stable  flow  aloft. 

b.  The  Mixed  Laver 

Deardorff  and  Willis  (1975)  predicted  diffusion  phenomena  in  convective  boundary 
layers  based  on  laboratory  experiments  and  sub-grid  scale  modeling  calculations. 
They  found  that  plume  behavior  was  completely  unlike  Gaussian  dispersion  models. 
Subsequent  modeling  approaches  for  dispersion  in  the  mixed  layer  are  reviewed  by 
Lamb  in  Nieustadt  and  van  Dop  (1984).  The  most  successful  models  have  adopted 
a  Lagrangian  viewpoint  and  use  sub-grid  scale  modeling.  Some  investigators  argue 
that  popular  approaches  of  Reynolds  stress  closure  do  not  work  well,  because  such 
models  fail  to  take  into  account  the  three-dimensional  effects  of  general 
subsidence  which  accompanies  the  rise  of  thermals  from  the  ground. 

As  stated  by  Lamb  (1984)  "the  long  life  times  of  updrafts  and  downdrafts  make  it 
very  frequently  possible  for  diffusing  particles  to  traverse  the  entire  depth  of 
the  mixed  layer  in  one  steady  motion.  This  is  incompatible  with  the  tenets  of 
Eulerian  K-theory,  in  which  turbulent  diffusion  is  envisaged  as  a  brownian  motion 
or  random  walk  type  process  in  which  particles  can  traverse  large  distances  after 
suffering  many  upward  and  downward  displacements.  .  for  this  reason,  dispersion 
parameters  derived  for  convective  conditions  from  K-theory  principles  are 
unreli able.  .  .Even  second  order  closure  models  are  likely  to  yield  erroneous 
predictions  of  concentration  in  the  mixed  layer."  Nonetheless,  Rodi  (1986)  has 
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recently  argued  kinetic-energy-closure  models  can  still  "effectively”  reproduce 
temperature  distribution  and  mixing-layer  depth  variations  with  time.  These 
features,  however,  will  not  be  sensitive  to  the  affect  of  subsidence  and  updrafts 
on  isolated  polluted  plumes. 

c.  The  Near-wall  Sublayers 

The  earliest  results  concerning  the  turbulence  structure  of  the  regions  beneath 
the  "mixed  layer"  are  due  to  Prandtl  (1932).  Later  other  researchers  such  as 
Obukhov  (1946,  1960),  Monin  i  Obukhov  (1954)  and  Priestly  (1954,  1960)  developed 
his  ideas  further  and  made  them  more  precise.  The  arguments  used  to  describe  the 
near-ground  flow  were  based  on  similarity  and  conventional  dimensional  arguments. 
They  suggested  that  shear,  temperature  and  turbulence  characteristics  could  be 
described  in  certain  layers  by  various  power-law  expressions  made  up  of  scaling 
variables  and  z  within  the  limits  of  unidentified  universal  constants.  Two 
layers  were  identified--a  near  wall  region  where  fluid  properties  were  expected 
to  vary  as  f,  =  z/L  when  (0  <  f  <  1): 

dU/dz  =  u,/(/cz)  0,(0,  dT/dz  =  -T./z  0,(f), 

=  U.^3(f), 

*  T,  04(f),  etc. 

and  a  free  convection  region  where  thermal  gradients  and  moments  were  expected 
to  be  independent  of  u*  when  f  »  1: 

dU/dz  =  8,u,2(q^)-i/3 2-4/3  ^  3^^-i/3^ 

dT/dz  =  or  0,  = 

a,  =  C,  (q^z)’'^  etc. 

Initially,  it  appeared  that  such  expressions  when  combined  with  mixed-layer 
scaling  described  the  unstably  stratified  environment  quite  well,  but  then 
several  inconsistencies  became  apparent: 

1.  The  free-convection  expressions  appeared  to  correlate  data  at 
unexpectedly  small  values  of  f  down  to  0.1, 

2.  Agreement  of  shear  and  temperature  gradients  with  the  derived 
expressions  were  rather  poor  as  f  got  larger  than  2  or  3, 

3.  It  was  not  clear  why  the  expression  for  dU/dz  in  the  free-convection 
layer  should  contain  u. ,  while  the  other  expressions  were 
independent  of  this  variable,  and,  finally, 

4.  Measurements  for  <ut>  and  <uwt>  did  not  agree  with  the  proposed 
expressions  at  all! 

Based  on  the  concepts  of  "directional"  or  "vector"  or  "anisotropic"  dimensional 
analysis  Betchov  and  Yaglom  (1971)  proposed  a  three-layer  model  for  the 
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atmospheric  surface  layer  which  lies  beneath  the  mixed  layer.  This  ''ess 
conventional  analysis  technique  proposes  that  horizontal  and  vertical  scales  of 
motion  and  velocity  should  be  scaled  by  separate  variables.  Therefore  it  permits 
one  to  obtain  sharper  results  than  those  implied  by  conventional  dimensional 
analysis.  In  some  cases  it  requires  an  increase  in  the  list  of  parameters  of  a 
problem  (the  additional  parameter  is  not  necessary  if  horizontal  and  vertical 
motions  are  energetically  uncoupled.) 

The  three  sublayers  and  their  limits  were  described  above.  Dimensional  arguments 
show  that  each  layer  requires  its  own  unique  shear,  temperature  and  turbulence 
parameterization  and  correlations.  In  the  dynamic  sublayer  where  buoyancy 
effects  can  be  neglected  the  well-known  logarithmic  formulae  for  U  and  T  are 
val id : 


U(z)  =  u,//c  ln{z/Zu),  T(z)  =  T.A  InCz/z^), 
u, ,  =  AjU,,  =  AjU, , 

Oj  =  A4T, ,  etc. 

where  T*  is  the  wall  temperature  and  z^  and  z^  are  surface  roughness 
coefficients. 

But  in  the  dvnamic-convecti ve  sublayer  the  velocity  scales  are  described  by  a 
vertical  scale  w,  ^  (Q^z) and  a  horizontal  scale  u,,  =  u,^/w, ,  and  a 
temperature  scale  0,  =  Q/w, ,  which  result  in  expressions: 

dU/dz  =  B,jU,./z  =  ByU.^  (Q;9)’'’^^z''‘''^  or 

dT/dz  =  -B,9./z  =  or  0,  = 

=  BjuJ  (Q^z)*’'^  a,  =  B3(Q^z)''^ 

a,  =  B,9.  =  B,Q2'^()9z)■^'^  etc. 

Then  in  the  free-convection  sublayer  only  w.  and  9,  are  required  such  that: 
dU/dz  =  C,w./z  =  or  0^  = 

dT/dz  =  -Ct9./z  =  or  0,  = 

=  C,  (Q^z)’'^  a,  =  C2(Q^z)’'^ 
etc. 

Unfortunately,  the  data  available  in  1971  was  insufficient  in  scope  or  quality 
to  confirm  these  expressions.  As  a  result  the  concepts  lay  dormant  until  1990 
when  Kader  and  Yaglom  reported  the  results  from  new  atmospheric  measurements  made 
in  Russia  during  the  mid-1980s  As  noted  in  Figure  1  derived  from  their  paper 
the  new  data  confirmed  the  non-monotonicity  of  the  function  0^,  which  had  never 
been  noticed  before.  They  uncovered  strong  evidence  for  the  existence  of  all 
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three  sublayers  and  even  proposed  values  for  many  of  the  universal  constants: 
A^j  =  2.5  ,  =  2.4,  By  =  1.7  ,  8^  =  1.1, 

Cjj  =  0.7  ,  C,  =  0.9, 

Ai  =  2.7  ,  Aj  =  2.5,  A3  =  1.25,  A^  =  2.9, 

83  =  1.65,  84  =  1.4 

C3  =  1.3  ,  C4  =  1.5. 


Indeed  Kader  and  Yaglom  were  even  able  to  suggest  c 
the  higher  third-order  moments  like  <w^t>,  <uvr>,  < 
rate,  e. 


coefficient  values  for  some  of 
uwt>,  etc.  and  the  dissipation 


The  behavior  of  the  functions  in  the  transition  regions  between  the  various 
sublayers  can  not  be  found  from  dimensional  reasoning,  but  empirical  formulae  can 
be  fitted  to  the  data,  for  example: 


2.5 


1  +  O.lf^  ■ 

.1  +  3  r  J  ’ 


0T(r) 


1.6 


3  +  r 

1  +  4f  +8r^ 


Such  expressions  can  be  integrated  over  a  range  of  f  values  to  produce  predictive 
expressions  for  the  velocity  and  temperature  profiles  as  noted  in  Figures  2  and 
3. 


d.  Entrainment  Interfacial  Laver 

The  rate  of  propagation  of. the  mixed  layer  into  the  stably  stratified  layer  above 
is  perhaps  the  least  well  understood  of  the  mixing  phenomena  which  occur  in  the 
CBL.  Since  the  process  takes  place  in  a  relatively  thin  layer  involving  the 
penetration,  propagation,  breaking  and  dissipation  of  internal  waves  it  has  been 
singularly  unavailable  to  accurate  measurement.  Nonetheless,  many  analytic  and 
numerical  models  have  been  proposed  based  on  simplistic  laboratory  measurements 
to  explain  the  phenomena.  Since  the  time  scales  of  this  process  are  rather  long 
comp'red  to  those  of  the  layers  below,  it  is  often  possible  to  study  the  two 
regions  separately.  Meteorologists,  oceanographers,  limnologists  and  solar  pond 
engineers  have  all  focused  on  interfacial  layer  behavior  through  such  simulations 
(Wyngaard,  1988;  Atkinson  and  Harleman,  1983;  Iwasa  et  al . ,  1988). 

Nonetheless,  such  assumptions  are  by  no  means  always  permitted.  Entrainment- 
induced  turbulence  can  extend  down  through  the  mixed  layer  to  the  surface, 
particularly  in  the  case  of  humidity  (Wyngaard,  1988).  The  entrainment  flux  at 
the  CBL  top  can  exceed  that  at  the  surface.  Thus  one  must  be  careful  when 
inferring  fluxes  of  momentum  or  heat  from  measurements  of  other  variables, 
especially  in  the  surface  layers. 

Most  laboratory  experiments  to  date  have  been  performed  in  small  apparatus  sized 
from  centimeters  to  a  meter  in  depth.  These  experiments  have  contributed  a  great 
deal,  but  they  leave  unresolved  questions  about  Reynolds  number  similarity 


5 


constraints  and  the  details  of  mixing  processes  too  small  to  resolve  at  such 
scales. 


2.  Laboratory  Simulation  of  the  CBL 

Much  of  the  understanding  of  the  CBL  has  come  from  convection  tank  and  large  eddy 
simulation  (LES)  studies.  The  original  work  by  Willis  and  Deardorff  (1974)  is 
now  complimented  by  wind-tunnel  measurements  taken  at  Colorado  State  University 
(Arya  and  Plate,  1969;  Poreh  and  Cermak,  1984),  and  the  soon  to  be  published  work 
by  Sawford  and  Hibbard  of  the  Atmospheric  Research  Division,  CSIRO  Aspendale, 
Australia.  Sawford  and  Hibbard  are  using  a  salt-water  analogue  to  the  thermal 
CBL  process  to  simulate  mixing  layer  phenomena. 

Thermally  stratified  wind  tunnels  currently  exist  at  Colorado  State  University, 
USA  (2  m  X  2  m  cross-section,  30  m  long)  and  at  the  National  Institute  for 
Environmental  Studies,  Japan  (  2  m  x  3  m  cross-section,  24  m  long).  New  wind- 
tunnel  facilities  specifically  designed  to  study  the  CBL  are  under  construction 
or  evaluation  at  the  Central  Electricity  Generating  Facility,  Leatherhead,  UK 
(1.5  X  3.5  m  cross-section,  20  m  long),  the  Ecole  Centrale  de  Lyon,  France  (2  x 
3.7  m  cross-section,  15  m  long),  Karlsruhe  University,  BRD  (0.5  x  1  m  cross- 
section,  4  m  long),  and  the  very  ambitious  facility  at  Monash  University, 
Australia  (5  m  x  10  m  cross-section,  40  m  long).  At  this  time  no  results  have 
been  published  based  on  experiments  performed  in  these  new  facilities. 

a.  Water-Tank  Experiments 

Laboratory  experiments  in  water  tanks  have  contributed  substantially  to  a  better 
understanding  of  convectively  mixed  layers.  Experiments  with  such  facilities 
continue  to  contribute  important  information  about  the  mixed  layer  and  its  rates 
of  growth.  Visualization  experiments  performed  with  laser  sheets  and  chemically 
fluorescing  plumes  yields  detailed  information  about  the  mixing  characteristics 
of  buoyant  plumes  in  the  CBL.  These  studies  clearly  demonstrate  that  the  process 
is  non-Gaussian. 

Small  size  limits  the  spatial  resolution  characteristics  of  most  water  tank 
simulation  of  the  CBL.  Typical  mixed  layer  depths  examined  are  from  0.1  to  0.5 
m  deep;  hence,  the  equivalent  atmospheric  surface  layer  regions  are  below  1  cm, 
and  the  stability  length,  L,  may  be  at  most  of  the  order  of  millimeters.  It  is 
unlikely  then  that  CBL  water  tank  measurements  will  contribute  substantially  to 
a  better  understanding  of  the  atmospheric  surface  layer  (ASL)  where  the  three- 
sublayers  reside. 

Deeper  entrainment  layers  may  be  produced  in  larger  water  tank  facilities.  For 
example  the  Stratified  Mixing  Facility  at  Colorado  State  University  consists  of 
a  2  m  deep,  10  m  long  tank  provided  with  thick  side  wall  insulation,  floor 
heating,  recirculation  pumps,  a  salt  water  fill ing  system,  visual  ization  windows, 
and  temperature  and  concentration  measurement  instrumentation. 
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b. 


Wind-Tunnel  Experiments 


Deardorff  (1972)  showed  for  the  CBL  that  when  -z,/L^  >  10  (Note:  Arya  (1982) 
suggests  >  5),  then  turbulence  scales  with  the  mixing  layer  depth,  ,  and  the 
convective  velocity,  w, .  Characteristic  distributions  of  temperature,  T, 
velocity,  U,  and  heat  flux,  <wt>,  are  sketched  in  Figure  4.  Time  is  expected  to 
scale  as  tw./z^  ,  vertical  distance  should  scale  as  z/z, ,  horizontal  distance  as 
xw,/(z^U^),  and  concentration  as  Cz^^U^/Q^  o>"  C,  z^  U^/Q^ .  Assuming  undistorted 
vertical  and  horizontal  scales,  this  means  one  must  maintain  model  and  prototype 
w,/U^  equivalent. 

Once  the  mixing  layer  depth  and  convective  velocity  scales  are  stipulated  for  the 
laboratory  experiment,  it  is  also  necessary  that  the  rate  of  mixing  layer  growth 
(dz^/dt)/w*  be  similar.  Unfortunately,  for  convenient  model  velocities,  eg., 

=  1  to  2  m/s,  very  large  surface  heat  fluxes  may  be  required,  eg.,  200  to  15.000 
watts/m^,  and  intense  stable  temperature  gradients  at  the  top  of  the  mixing 
layer,  eg.  dT/az  =  50  to  2000  *C/m.  An  alternative  approach  to  the  creation  of 
an  elevated  inversion  which  might  be  adequate  is  to  use  the  upper  wind-tunnel 
roof  as  the  effective  inversion  height,  z^  .  The  validity  of  this  approach  must 
be  examined  carefully,  however,  since  a  downward  flux  of  heat  due  to  interfacial 
entrainment  plays  an  important  role  in  many  mixed  layer  characteristics. 


3.  Laboratory  Simulation  of  Atmospheric  Transport  in  tne  CBL 

Simulation  of  the  behavior  of  industrial  plumes  released  into  the  CBL  require  the 
additional  equality  of  plume  exit  momentum  ratio, 

MR  =  A3W32  d,V(pU'z,2), 

and  the  ratio  of  plume  exit  buoyancy  to  frea-stream  momentum, 

Fr^  =  gApD3W3/(pU^ ) , 


sometimes  called  the  modified  Froude  number.  The  joint  stipulation  that  the 
velocity  ratio,  w,/U,  length  ratio,  D^/Zj  ,  and  modified  Froude  nu,  ^er,  Fr^,  are 
similar  is  equivalent  to  equality  of  the  buoyancy  parameter,  F. ,  defined  as  the 
ratio  of  plume  exit  buoyancy  to  free  stream  convective  buoyancy. 


z^  Uw,^ 

and  for 


4p 


n\ 

Ds  1 

■  U  ■ 

4 
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F.  was  first  used  to  describe  plume  motion  in  C8L  situations  by  Briggs  (1975). 
but  it  has  also  been  affirmed  by  Kerman  (1979).  Willis  and  Deardorff  (1983)  and 
Venkatram  (1991).  When  F,  is  zero,  plume  trajectory  and  dispersion  is  dominated 
by  convective  turbulence,  at  0.02  weak  buoyancy  interactions  occur,  at  0.05 
moderate  buoyancy  effects  are  observed,  and  by  0.10  strong  buoyancy  effects  can 
be  expected. 

An  attractive  characteristic  of  the  buoyancy  parameters  Fr^  or  F,  is  that  they 
permit  the  enhancement  of  model  velocities  through  distortion  of  plume  density 
Indeed,  using  pure  helium  as  a  model  plume  simulant  allows  one  to  increase 
typical  model  flow  velocities  2.5  times.  Model  plume  Reynolds  number  also 
increases  by  a  factor  of  two,  which  encourages  the  emission  of  a  turbulent  model 
plume.  The  only  distortion  of  scale  caused  by  this  modeling  is  in  the  kinematic 
ratios  near  the  plume  exit,  and  mass  ratios,  which  need  consideration  when 
determining  concentration  ratios  between  model  and  full  scale. 

a.  Exhaust  characteristics  of  typical  power-station  plumes 

Two  typical  prototype  power  stations  are  considered:  the  Gladstone  Power  Station 
and  the  Loy  Yang  or  Driffield  Power  Station,  Australia.  Consider  a  sunny  day 
with  an  inversion  height  of  1000  m,  wind  speeds  of  U  =  4  ms  \  and  an  ambient 
temperature  of  T  =  20*0.  Take  Q  =  400  Wm'^ ,  which  corresponds  approximately  to 
Pasquill  stability  class  B.  Their  characteristics  are  summarized  in  Table  1. 
Based  on  calculated  F,  values  the  Gladstone  Station  plume  will  be  strongly 
affected  by  convective  turbulence,  whereas  the  Loy  Yang  Station  plume  will  be 
dominated  initially  by  plume  buoyancy. 

Table  1  Typical  Power-station  Exhaust  Characteristics 


Power  station: 

Gladstone 

Loy  Yang 

Stack  height,  h^  (m)  = 

140 

260 

Stack  diameter,  (m)  = 

6 

11 

Emission  velocity,  (m  s'M 

15 

23 

Emission  temperature,  T^  (*C) 

140 

175 

Emission  density,  (kg  m‘^ ) 

0  851 

0.785 

Stack  density  defec*,  Ap  (kg  m'^ ) 

0.349 

0  415 

Convective  velocity,  w,  (m  s  ’) 

2,23 

2.23 

Buoyancy  flux,  F,  (m‘^  s'^) 

385 

2874 

Buoyancy  parameter ,  F. 

0.0193 

■H 
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b. 


Consider  a  wind-tunnel  simulation  where  a  helium  simulant  gas  is  specified, 
maximum  model  surface  heating  rate  is  5  KWm'^ ,  and  three  model  scales  750:1, 
500:1  and  200:1  are  possible.  Assume  one  examines  behavior  10  km  downwind  of  the 
stack  for  the  atmospheric  conditions  given  in  section  3-a  above.  In  such 
circumstances  the  similitude  parameters  discussed  above  reduce  to  the  model  data 
in  Table  2. 


Table  2  Model  conditions  to  simulate  Gladstone  and  Loy  Yang  Power  Station 
exhaust  plumes. 


Model  Parameter 

Gladstone  Model 

Loy  Yang  Model 

750 

500 

200 

750 

500 

200 

X  (m)  = 

13.3 

20 

50 

13.3 

20 

50 

z,  (m)  = 

1.33 

2 

5 

1.33 

2 

5 

h=  (m)  = 

0.19 

0.28 

0.70 

0.35 

0.52 

1.30 

D„  (m)  = 

0.008 

0.012 

0.030 

0.015 

0.022 

0.055 

0.38 

0.46 

0.73 

0.34 

0.42 

0.66 

w.  (ms'M  = 

3.21 

3.94 

6.22 

5.18 

6.34 

10.05 

q  (Wm'’  )  = 

254 

312  J 

490 

184 

225 

356 

mm^m 

0.21 

0.26 

0.41 

0.19 

0.23 

0.37 

w./u  = 

8.49 

8.49 

8.49 

15.24 

’  15.24 

15.24 

w./U  = 

0.56 

0.56 

0.56 

0.56 

0.56 

0.56 

F..  (mS-^)  = 

.00043 

0.0012 

0.0118 

0.0024 

0.0065 

0.064 

Fr™  = 

4.01 

4.01 

4.01 

16.3 

16.3 

F.  = 

0.019 

0.019 

0.019 

0.145 

0.145 

MR  = 

.00036 

.00036 

.00036 

.00039 

0.0039 

0.0039 

Re,  = 

214 

394 

1555 

646 

1162 

4610 

« 

X  = 

5.58 

5.58 

5.58 

5.58 

5.58 

5.58 

If  one  requires  model  mixed  layer  velocities  greater  than  0.5  ms’’  and  limits 
the  wind  tunnel  length  to  40  m;  then  wind  tunnel  models  of  both  the  Gladstone  and 
the  Loy  Yang  situations  are  practical  at  length-scale-ratios  between  about  400:1 
to  300:1.  As  required  the  model  and  prototype  values  of  convective  veloci  ty 
ratio  (w./U),  stack  hei^ght  ratio  (hj/z, )  ,  buoyancy  parameter  (F,),  and 
dimensionless  distance  (x*  =  xw./(z,U))  will  all  be  equivalent. 
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4.  Wind-tunnel  Operating  Range  for  C8L  Simulations 

Meteorological  wind  tunnels  are,  in  effect,  analog  computers  with  'near¬ 
infinitesimal'  resolution  and  'near-infinite'  memory  (Snyder.  1972).  They  employ 
real  fluids,  not  mathematical  models  of  fluids,  and  produce  inherently  viscous, 
turbulent,  nonhydrostatic,  non-Boussinesq,  and  compressible  flows  with  no-slip 
boundary  conditions.  However,  flows  in  scaled  physical  models  are  also  only 
partially  similar  and  cannot  at  present  include  all  processes  present  in  the 
atmosphere  such  as  Coriolis  acceleration,  exchange  of  energy  by  radiation, 
conduction  into  the  soil,  and  phase  changes  of  water. 

Simulation  of  atmospheric  motions  by  wind-tunnel  flows  has  occurred  for  almost 
100  years  since  Professor  LeCour  constructed  a  wind-mill  test  facility  in  Askov, 
Denmark,  in  1895  and  Gustaf  Eifel  designed  his  exhibition  tower  in  1889. 
Background  reviews  about  laboratory  simulation  were  prepared  by  Cermak  (1975), 
Davenport  and  Isyumov  (1967)  and  Melbourne  (1977).  Meroney  (1981,  1990) 
considered  the  simulation  of  complex  terrain  and  valley  drainage  situations. 
Snyder  (1981)  suggested  similarity  criteria  for  the  study  of  air-pollution 
meteorology  in  near  neutral  situations.  Meroney  (1987)  extended  the  discussion 
to  the  simulation  of  dense-gas  plumes  in  the  surface  layer.  Meroney  et  al . 
(1975)  and  Avissar  et  al .  (1990)  proposed  simulation  criteria  and  operating 
ranges  for  the  simulation  of  sea  and  land  breezes. 

a.  General  similarity  requirements 

During  physical  model  simulations  of  atmospheric  flows,  scale-model  replicas  of 
observed  ground-level  buildings  and  terrain  are  constructed  and  inserted  into  a 
laboratory  flow  facility.  The  flow  characteristics  and  stratifications  of  the 
air  in  the  wind  tunnel  are  adjusted  to  be  similar  as  possible  to  the  atmospheric 
conditions.  Complete  equivalence  of  the  laboratory  model  and  atmospheric 
prototype  flow  fields  requires  geometric,  kinematic,  dynamic,  and  thermal 
similarity.  In  addition,  boundary  conditions  upstream,  downstream,  at  the  lower 
surface,  and  near  the  top  of  the  physical  model  must  be  similar  to  those  at  the 
corresponding  boundaries  of  the  modeled  atmospheric  domain.  These  multiple 
similarity  requirements,  the  characteristics  of  the  wind  tunnel  and  its 
instrumentation,  and  the  nature  of  the  atmospheric  phenomenon  to  be  modeled  all 
help  to  determine  the  operating  range  (OR)  for  a  wind-tunnel  simulation. 

Similarity  characterization  of  stratified  atmospheric  flows  are  summarized  by 
Avissar  et  al .  (1990).  Specific  characterization  of  CBL  phenomena  was  discussed 
earlier  in  Sections  1.1  through  1.4.  Equality  of  these  similitude  parameters 
must  be  supplemented  by  the  requirements  that  the  surface  boundary  conditions  and 
the  approach- flow  characteristics  also  be  similar  for  model  and  prototype. 
Boundary-condition  similarity  requires  similar  values  of 

0  Surface  roughness, 

0  Topographic  rel ief , 

0  Surface  temperature  distribution, 

0  Upstream  distribution  of  mean  and  turbulent  velocities, 

0  Upstream  distribution  of  mean  and  turbulent  temperatures,  including 
inversion  height,  and 
0  Longitudinal  pressure  gradient. 
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If  all  of  the  above  conditions  are  met  simultaneously,  then  all  scales  of  motion 
ranging  from  the  atmospheric  microscale  to  mesoscale  could  be  simulated  exactly 
by  the  laboratory  model.  Unfortunately,  not  all  conditions  can  be  satisfied 
simultaneously  by  a  scaled  model  since  some  are  incompatible  or  conflicting; 
hence,  only  partial  or  approximate  similarity  can  be  achieved.  This  suggests 
that  a  laboratory  model  for  a  particular  meteorological  situation  must  be 
designed  to  simulate  most  accurately  those  scales  of  motion  which  are  of  greatest 
significance  for  the  application.  In  the  case  of  a  CBL  with  undistorted 
horizontal  and  vertical  scales  the  primary  parameters  of  interest  will  be  w./U 
and  inversion  height,  z^  . 


b.  Instrumentation  Characteristics 

Avissar  et  al .  (1990)  reviewed  wind-tunnel  instrumentation  characteristics,  and 
they  suggested  that  measurement  accuracy  was  enhanced  for  model  flow  speeds  of 
0.5  m/s  and  greater.  Averaging  times  and  sampling  rates  may  be  expected  to 
produce  fractional  errors  of  1,  5  and  10%  for  averaging  times  of  400,  16,  and  4 
seconds,  respectively.  Conventional  sized  hot-film  instrumentation  should 
produce  only  a  4%  spatial  resolution  error. 

c.  Wind-tunnel  characteristics 

Simulation  of  the  atmospheric  CBL  is  not  only  a  function  of  the  governing  flow 
physics  but  also  depends  on  the  availability  of  a  suitable  simulation  facility 
and  its  instrumentation.  In  particular  the  size  of  the  wind-tunnel  test  section 
will  determine  the  smallest  model  length  scale  ratio  (LSR).  Most  meteorological 
tunnels  range  in  size  from  0.5  x  0.5  m  to  those  with  working  cross-sections  of 
3  X  4  m.  An  exception  is  the  new  environmental  wind  tunnel  at  Monash  with  a 
cross-section  of  5  x  10  m.  Density  stratification  can  be  induced  by  use  of  heat 
exchangers,  injection  of  heated  air,  gases  of  different  molecular  weight,  or 
latent  heat  absorption  or  release  during  phase  change  (e.g.,  Ogawa  et.  a/.,  1985; 
Meroney,  1986).  By  using  vortex  generators,  fences,  roughness  elements,  grids, 
screens  or  jets,  a  wide  range  of  turbulence  integral  scales  can  be  introduced 
into  the  tunnel  boundary  layer.  Choice  of  model  surface  roughness  or 
stratification  permits  control  of  surface  turbulence  intensity,  dimensionless 
wall  shear,  and  velocity  profile  shape. 

Two  tunnels  have  been  chosen  to  represent  the  characteristics  of  such  facilities. 
The  meteorological  wind  tunnel  (MWT)  at  Colorado  State  University  is  a  large, 
closed-circuit  facility  with  a  1.8  m  high  by  1.8  m  wide  by  24  m  long  test 
section.  Wind  speeds  are  continuously  variable  from  0.1  to  30  m  s'’  and  ambient 
air  temperatures  can  be  varied  from  5  to  205’C  (e.g.  Cermak,  1975,  1982).  Ten 
meters  of  the  upstream  floor  can  be  cooled  between  I’C  and  ambient  temperature 
while  12  m  of  the  downstream  test-section  floor  can  be  heated  from  1  to  200’C. 

The  environmental  wind  tunnel  (EWT)  at  Monash  University  is  a  large,  open-circuit 
facility  with  a  5-7  m  high  by  10-12  m  wide  by  40  m  long  test  section.  Wind 
speeds  are  continuously  variable  from  0.1  to  18  ms’’  (  or  to  45  ms'’  with  cross- 
section  reduced)  .  Ambient  air  temperatures  are  drawn  into  the  entrance  section, 
but  wall  heaters  may  be  inserted  to  regulate  wall  temperatures.  The  wall  heaters 
have  a  surface  heat  capacity  of  5  kw  m'^ . 


II 


The  following  criteria  have  been  identified  which  limit  the  MWT  and  EWT 
operational  ranges: 


Colorado  State  Monash 
_ _ EWT 

(?  Maximum  model  inversion  height  <  1  m  <  5  m 

(?  Maximum  inversion  strength  <  -10*Cm‘’  <  -10*Cm'’ 

(?  Minimum  model  inversion  height  0  0 

@  Maximum  model  blockage  <  5  %  <  5  % 

(?  Minimum  plume  Reynolds  number  >  300  >  300 

@  Minimum  model  wind  speed  >  0.5  'ns’’  >  0.5  ms’’ 

@  Maximum  model  heat  flux  <5  KWm'^  <  5  KWm’^ 


Coriolis  force  considerations  also  limit  the  maximum  acceptable  LSR.  Snyder  (1972) 
suggests  a  5  km  maximum  cut-off  point  for  horizontal  length  scales  for  modeling 
atmospheric  diffusion.  Mery  et.  al .  (1974)  suggests  a  15  km  limit,  Ukeguch:  et  al . 
(1967)  suggest  a  40  to  50  km  limit,  and  Cermak  et  a/.  (1966)  proposed  a  150  km 
limit.  Given  the  strong  mixing  present  in  CBL  situations  and  time  scales  of  about 
3  hours  then  surface  generated  stresses  should  dominate  most  situations  for  at  least 
2  to  5  km. 


The  most  severe  restrictions  on  the  wind-tunnel  operating  range  or  performance 
envelope  result  from  geometric  similarity  constraints.  For  example  Deardorff 
proposed  that  laboratory  conditions  should  be  sought  where  z,/L  >  10  (z^L  >  25) 
to  simulate  the  mixed  layer.  Yet  the  height  of  most  facilities  (1-2  m)  would  limit 
the  dynamic  layer  to  2  mm,  the  dynamic  convective  layer  to  2  cm,  and  the  free 
convection  layer  to  less  than  20  cm.  Presuming  wind-tunnel  facilities  with  model 
z,-  =•  1  and  5  m,  respectively,  then  the  various  characteristic  model  depths  might 
be: 


Colorado  State  Monash 
MWT  _ EWT 


~  1  m  5m, 

U  =  0.002  m  0.01  m, 

L  =  0.02  m  0.10  m, 

L,,  =  0.04-0.10  m  0.20-0.50  m,  and 

=  0.04-0.12  m  0.20-0.60  m. 

An  additional  criteria  required  to  avoid  secondary  circulations  imposed  by  facility 
size  is  that  the  aspect  ratio  of  width  to  inversion  height  should  be  greater  than 
4.  Willis  and  Deardorff  (1974)  found  that  at  values  of  2,  results  were  no  longer 
consistent  with  an  infinitely  large  homogenous  layer. 

Based  on  the  previous  discussions,  the  following  seven  similarity  criteria  appear 
pertinent  to  the  physical  simulation  of  stack  plumes  dispersing  within  CBL 
circulations: 
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*  [R^bU  above  the  inversion; 

2.  [w*/U]^  =  [w*/U]p  in  the  mixed  layer; 

3.  [h,/zj^  =  [h3/z,]p; 

4.  [PsCWjVCpU')],  =  [P3(Ws)V{pU2)]p; 

5.  [gApD3W^/(pU^)]„  -  [gA/)0,W,/(pU^)]p; 

6.  (Re)  piuw  ^  ^^fnin  300  ^  and 

7.  Similar  upwind  velocity,  temperature,  and  turbulence  profiles. 

In  a  CBL  situation  mechanical  and  thermal  turbulence  influence  one  another  in  the 
surface  layers.  In  turn  this  alters  the  velocity  and  temperature  profiles 
throughout  the  CBL.  The  bulk  Richardson  number  equality  (criterion  1)  requires  that 
above  the  CBL  inversion 


[(5T/az),  1 

1/2 

1 

1/2 

1 

.{5e/3z)p  , 

T, 

where  LSR  denotes  the  prototype-to-model  vertical  length-scale  ratio  and  subscripts 
m  and  p  refer  to  model  and  prototype  conditions,  respectively. 


The  requirement  that  the  convective  velocity  to  advection  velocity  scale  ratio 
(Criterion  2)  is  also  satisfied  implies  that 


LSR 


(w.)p 


Similarity  of  the  depth  of  the  CBL  (Criterion  3)  requires  that 


LSR  =  "  ''P  . 

(Zi). 

Stipulation  that  the  modeled  plume  momentum  and  buoyancy  conditions  (Criterions  4 
and  5)  are  similar  results  in  the  requirement  that 


LSR 


[(Ps  -  ' 

1/2 

[“p  1 

iPs  -  P)p  . 

ip.p)^  . 

The  plume  Reynolds  number  will  be  large  enough  to  ensure  a  fully  turbulent  plume 
(Criterion  6)  if  one  stipulates  model  conditions  such  that 


(Ws)o(Ds)o 
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Combining  the  similarity  constraints  given  by  the  equations  for  LSR  above  with  the 
characteristic  sizes  and  flow  capacities  of  the  wind  tunnels  described  in  Section 
4-c  and  typical  atmospheric  CBL  and  power-plant  exhaust  conditions  noted  in  Section 
3  provides  the  relationships  and  data  needed  to  construct  the  wind-tunnel  operating 
range  and  to  identify  reasonable  simulation  scenarios.  Figures  5  through  10  aisplay 
performance  envelopes  for  the  two  wind  tunnel  sizes  that  are  required  to  produce 
Pasquill  category  A,  B,  and  C  style  CBL  stratification,  respectively.  The  shaded 
areas  are  excluded  per  criteria  limitations  for  both  CBL  and  plume  rise  simulations. 
Model  plumes  produced  under  CBL  conditions  in  the  cross-hatched  region  may  not  be 
turbulent.  Simulation  of  the  Gladstone  or  Loy  Yang  power  station  plumes  must  lie 
along  the  lines  and  within  the  remaining  open  polygon  region. 

Although  the  smaller  MWT  can  marginally  reproduce  CBL  conditions  (Poreh  and  Cermak. 
1984:  z,/L  =  5.6-10.7),  simultaneous  simulation  of  plume  behavior  as  stipulated 
by  MR  and  F,  will  be  unlikely  if  not  impossible.  Even  for  the  EWT  joint  simulation 
of  the  CBL  and  stack-plume  behavior  can  only  be  performed  over  a  limited  LSR  and  U.^ 
range,  especially  for  Pasqui 1 1 -Gi f ford  A  stability  situations. 

It  is  concluded  that  an  extensive  range  of  wind-tunnel  modeling  of  the  dispersion 
of  large  power  station  and  chemical  plant  plumes,  in  a  highly  unstable  atmospheric 
boundary  layer  could  be  achieved  with  a  wind-tunnel  working  section  5  m  high  by  10 
m  wide  by  40  m  long  with  1  MW  of  floor  heating.  Model  length  scale  ratios  would 
center  about  1/200,  with  the  smallest  being  about  1/400  and  the  largest  1/100.  The 
limitations  at  the  small  scale  are  Reynolds  number  and  the  practicability  of 
controlling  very  low  wind-tunnel  speeds.  The  limitations  at  the  large  scale  are 
mainly  the  enormous  heat  requirements  and  size  of  the  facility  needed  to  model 
adequate  downstream  distances  from  the  source. 


5.  Boundary- layer  Control  Methods  for  Augmenting  CBL  Simulation 

A  special  laboratory  facility  designed  to  produce  CBL  simulation  at  higher  wind 
speeds  would  be  very  desirable.  E.J  Plate,  Karlsruhe  University,  has  suggested 
that  it  may  be  possible  to  reduce  the  intense  surface  heat  flux  required  in  the 
laboratory  with  a  combination  of  wall  heating,  mass  transp ' rati  on  of  hot  fluid  at 
the  wall,  and  induced  subsidence  of  the  mixing  layer  by  side-wall  suction.  In  the 
following  paragraphs  we  will  examine  the  conditions  required  to  perform  augmented 
simulation  that  might  successfully  simulate  the  convective  boundary  layers  measured 
by  Kaimal  et  al .  (1976)  during  field  experiments  in  Minnesota,  U.S.A. 


a .  A  one-dimensional  model  for  augmented  CBL  simulation 

To  evaluate  the  viability  of  boundary  control  for  CBL  simulation  consider  rhe  one¬ 
dimensional  entrainment  models  first  proposed  by  Lilly  (1979)  and  later  improveo  by 
Betts  (1974)  and  Oeardorff  (1979).  They  proposed 

-<wt>,  =  (dz^  /dt  -  Wj,  )  A9 

where  the  heat  flux.  <wt>,  is  the  negative  heat  flux  of  entrainment  at  ?  Now 
Oeardorff  et  a/.  (1980)  has  shown  that  -<;wt>,  =  0  2<wt>^  .  Of'en  A0  is  of  the  order 
FAz,  and  Az,  is  of  the  order  of  0  2z, .  If  one  assumes  that  <wt>  =  <wt>  = 
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<wt>  +  Wy(T  -  Tj ) ,  where  is  a  wall  mass  transfer  rate  of  air  injected  at 
temperature  i then 

<wt>  +  W^{T^  -  TJ  =  (dZj/dt  -  W.jrZj,  or 
dz^/dt  =  [<wt>  +  W^(T^  -  TJ]/(rzJ  +  or 

(dz,/dt)/w,  =  0.3/(Ri.),,,  + 

where  (Ri.)eff  - 

Ri.  =  g^rz,.V[g^<wt>z,-]^^^, 

Ri,  =  g^rz,V[g^WJT,  -  TJ  z,]2/3^  and 

W.r  =  W,,/w.. 

In  order  to  maintain  the  correct  mixing  layer  growth  rate  in  the  model  one  must 
require 

[(dz,/dt)/w.]p,,,,,^  =  [(dz,/dt)M  - 

When  the  prototype  mixing  layer  maintains  a  constant  height,  then  W^.*  is  determined 
by  (Ri»)eff>  which  can  be  produced  by  various  combinations  of  the  parameters  Ri,  and 
Ri  .  Figure  11  displays  the  values  of  these  parameters  required  to  produce  various 
values  of  [{dz^/dt)/w,  Figure  12  specifies  combinations  of  these 

parameters  which  produce  zero  mixing  layer  growth  for  different  effective  Richardson 

numbers,  (Ri.)^^^. 

Unfortunately,  the  action  of  mass  addition  at  the  wall  which  helps  to  simulate  the 
mixing  layer  region,  also  affects  the  profiles  of  mean  velocity  and  temperature  in 
the  surface  layer,  *  O.lz,-.  As  noted  by  Arya  (1982),  the  surface  layer  is 

expected  to  be  a  function  of  both  surface  and  mixing  layer  scales.  The  surface 

layer  scales  are  the  friction  velocity,  u. ,  and  the  friction  temperature,  9,,  which 
are  related  to  the  mixing  layer  scales  by  the  expressions 


w./u.  -  and 

T./9.  - 

where  k  is  the  von  Karman  constant  equal  to  0.4.  Thus,  surface  layer 
cnaracteristics  are  expected  to  be  functions  of  all  of  these  variables. 

Mass  addition  to  unheated  boundary  layers  results  in  distortion  of  the  velocity  and 
temperature  profiles  near  the  wall  due  to  longitudinal  momentum  deficits.  Meroney 
(1967,  1968)  has  derived  expressions  for  the  velocity  and  temperature  profiles  in 
transpired  flow  which  obey  the  equations  of  motion  near  the  wall.  One  begins  with 
the  one-dimensional  equations  of  motion  with  mass  addition,  or 

W(dU/dz)  =*  d<-uw>/dz, 
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W(dW/dz)  =  d<w^>/dz  +  gAp/p, 
W(dT/dz)  =  d<-wt>/dz,  and 
dW/dz  =  0.0. 


The  continuity  equation  pguires  that  U  =  =  constant.  Prandfl's  mixing-length 

model  suggests  <-uw>  »  /c^z"  [dU/dz]  (dU/dz) .  One  can  extend  this  to  heat  flux  by 
<-wt>  =»  /c^z^  |dU/dz|  (dT/dz) .  Introducing  these  closure  models  into  the  equations 
above  and  integrating  once  produces, 

«;z(dU/dz)/u,  =«  (1  +  and 

/cz(dT/dz)/9.  *  (1  +  W"T")/(1  + 

where  U*  *  U/u, ,  W*  =  W  /u,  and  T*  =  {T  -  T  )/0, .  Meroney  (1967,  1968)  developed 
expanded  expressions  for  the  transpired  ooundary  layer  from  these  equations 
including  a  wake  correction  and  a  sublayer  term. 

In  the  near  wall  region  proposed  by  Monin-Obukhov  similarity  it  is  conventional  to 
express  the  dimensionless  shear  and  heat  flux  in  the  surface  layer  as 


Kz(dU/dz)/u.  »  (^ij(z/L^),  and  /cz(dT/dz)/e.  =  0^  (z/L^ ) , 

where  is  the  Monin-Obukhov  similarity  length.  For  unstable  flows  where  is 
negative  a  good  approximation  to  the  functionals  is 


0,  =  0/  =  (1  -  15(z/L^)) 


Combining  the  equations  above  to  account  for  the  joint  effects  of  surface  heating 
and  mass  transpiration  gives, 


/cz(dU/dz)/u.  =  (1  +  W*U*)''V(1  -  15{z/L  and 


Kz(dT/dz)/9. 


(1  +  wn*) 


(1  +  W*U*)’'2(i  -  15(z/L^))’''’ 

These  expressions  can  be  integrated  to  produce  the  following  functions, 
U/-U"  -  -WV(4k^)F^  +  (1  +  and 

T"  =  [(1  +  W*T/)exp(-G)  -  1]/W*, 


where 


F 

G 


I 


[x(l  -  15(L./L™)x)’'M-'  dx,  and 


-ASL  '  ‘-mo 


z/Lasl 


.V.  |j 


[x(l-15(L„,/L^)x)^^^((UW*U,*)’'^-v^*F/(2M] 

z/Lasl 


- 1 


dx  - 
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These  expressions  are  integrable  for  various  combinations  of  W* ,  (1^5^  ''-  )■  ' 
and  for  the  magnitudes  of  U*  and  T"  at  various  heights,  z,  ™ 

Figure  13  displays  the  variation  of  dimensionless  velocity,  U^*-U*,  versus  surface 
layer  position  for  finite  wall  blowing  and  surface  heating  rates.  Heating  tends  to 
increase  turbulence  and  reduce  the  velocity  profile  power-law  coefficient.  Figure 
14  displays  the  dimensionless  temperature  variation  for  finite  wall  blowing  and 
various  surface  heating  rates.  Blowing  tends  to  reduce  velocities  near  the  wall  and 
increase  the  velocity  profile  power-law  coefficient. 

b.  Augmented  wind-tunnel  conditions  to  simulate  the  CBL 

Results  from  a  CBL  experiment  conducted  over  a  flat  site  in  northwestern  Minnesota. 
USA,  are  described  by  Kaimal  et  al .  (1976).  Eleven  sets  of  data  are  tabulated  for 
CBL  conditions  found  during  the  afternoon.  For  one  set  of  data  the  values  of  , 

,  T3 ,  and  dz^dt  can  be  estimated  from  figures  provided.  For  the  other  cases 
assumed  values  for  these  parameters  must  be  used.  Table  3  summarizes  the  prototype 
conditions  considered.  To  calculate  missing  values  a  constant  surface  drag  of 
0.0014  was  assumed.  Missing  T^  and  T^  and  b  values  were  set  to  295°K,  300*K  and 
0.004  *Km/s,  respectively.  When  missing,  dz^/dt  was  set  to  zero.  The  similarity 
rules  discussed  above  permits  one  to  calculate  model  values  for  heat  flux,  mass  flux 
and  inversion  temperature  gradient,  T,  once  model  velocities,  temperatures  and 
inversion  height  are  specified. 

Use  of  the  mass  addition  option  can  reduce  the  wall  heating  considerably  for  only 
modest  wall  mass  transport.  If  one  assumes  the  maximum  convenient  wall  surface  flux 
is  200  w  m'^  ,  then  blowing  rates  of  only  2  to  20  cm  s’’  are  required  for  U^  =  2  m 
s’’  . 

The  above-mixing-layer  temperature  gradients  were  presumed  specified  by  the 
stipulation  that  (Ri.)eff  must  be  equal  to  the  prototype  value.  This  similitude 
requirement  results  in  unfortunately  large  temperature  gradients  for  the  higher 
model  wind  speeds.  Some  may  be  obtainable  in  the  laboratory,  others  are  definitely 
too  large.  But  the  values  required  for  subsidence  velocity,  W^^  ,  are  generally  less 
than  1  cm  s’’.  This  can  be  obtained  by  only  mild  suction  at  the  wall,  or  ignored 
entirely.  An  alternative  would  be  to  use  reduced  values  for  model  F,  but  definitely 
suck  at  the  walls  of  the  tunnel  to  reduce  the  expected  increase  in  mixing  layer 
growth . 

The  effect  of  mass  addition  on  the  simulated  surface  layer  profiles  may  indeed  be 
significant.  Bar  charts  which  indicate  the  order  of  magnitude  of  conditions 
required  to  simulate  the  Minnesota  experimental  cases  in  a  wind  tunnel  with  a  one- 
meter  height  are  shown  in  Figures  16  through  19.  A  comparison  of  profiles  for  Case 
2A1  are  found  in  Figure  18.  Near  the  wall  heating  will  reduce  the  effects  of  mass 
addition  to  a  nearly  neutral  profile  shape,  but  the  dimensionless  velocities  would 
still  be  larger  than  their  prototype  counterpart.  Temperature  profiles  of  T'-T^* 
display  a  similar  behavior  with  heating  and  blowing.  When  =  (R^)Lasi  =  0  then 
temperature  and  velocity  obey  the  same  equations  of  motion  and  boundary  conditions; 
hence,  U,*-U*  =  T„*-T*. 

’  rn 
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6.  Conclusions 


This  paper  has  examined  the  characteristics,  capabilities,  and  limitations  of 
meteorological  wind  tunnels  to  simulate  the  atmospheric  convective  boundary  layer. 

A  significant  conclusion  is  that  such  facilities  provide  an  opportunity  to  explore 
CBL  sublayer  and  entrainment  layer  behavior  and  an  opportunity  to  examine  routine 
interaction  of  plumes  with  CBL  turbulence. 

Extensive  experience  with  stratified  water  tanks  and  wind  tunnels  definitely 
suggests  that  the  important  turbulence  characteristics  of  the  convective  boundary 
layer  can  be  simulated  in  the  laboratory.  Nonetheless,  laboratory  simulation  is 
often  not  automatic  or  convenient.  Examination  of  the  characteristics  of  mixing 
layer  entrainment  and  surface  layer  behavior  determines  that: 

1.  Useful  CBL  simulations  may  be  obtained  in  sufficiently  large  stratified  wind- 
tunnel  facilities  without  augmentation  techniques. 

2.  The  use  of  sufficiently  large  stratified  wind-tunnel  facilities  will  provide 
a  means  to  study  the  atmospheric  sublayers  associated  with  the  CBL. 

3.  Large  to  modest  wind  tunnel  mixing  layer  velocities,  1  <  <  2  m  s’\  can  be 

obtained  in  smaller  facilities  through  the  use  of  floor  mass  addition  and 
wall  suction;  however, 

4.  Large  above-mixing-layer  model  temperature  gradients  may  be  required,  unless 

5.  Roof  mass  addition  and  side-wall  suction  is  used  to  reduce  the  model 
temperature  gradients  above  the  mixing  layer. 

6.  Mass  addition  and  boundary  layer  control  required  on  the  side  walls  are 
generally  modest  requiring  cross-flow  velocity  magnitudes  between  0.1  to  25 
cm  s'^  . 

7.  Mass  addition  will  modify  mean  velocity  and  temperature  profiles  in  the 
surface  layer.  The  additional  wall  heat  flux  due  to  mass  addition  is  not 
likely  to  completely  compensate  for  the  momentum  deficit  effects  of  mass 
addition. 

8.  The  complication  of  such  boundary-layer-control  structures  can  be  avoided  if 
one  utilizes  very  low  model  mixing  layer  velocities,  <  0.5  m  s’’.  An 
alternative  approach  would  be  to  develop  hot-film  technology,  laser  specal 
velocimetry,  or  1 aser-anemometry  to  make  such  measurements  in  the  lower 
velocity  range. 
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Table  3:  Prototype  wind  conditions  during  Minnesota  Convective  Boundary  Layer 
Experiments  of  Kaimal  et  ai . ,  1976 
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Figure  l  Dimensionless  velocity  and  temperature 
gradients  0,j  (f)  and  0^  (f )  obtained  from  Kader  and 
Yaglom  (1991) — lines  are  power  law  expressions. 
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Figure  2  Dimensionless  velocity  profiles  based  on 
Kader  and  Yaglom  (1991)  empirical  sublayer 
equations  for  dimensionless  shear. 
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Figure  3  Dimensionless  temperature  profiles  based 
on  Kaglom  and  Yader  (1990)  empirical  expressions 
for  dimensionless  temperature  gradient. 


Figure  4  Schematic  of  the  structure  of  a  typical 
convective  boundary  layer  (CBL) 
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Figure  7  Operating  range  for  CBL  simulation  of 
Pasquill  Category  C  conditions  in  EWT 
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Figure  8  Operating  range  for  CBL  simulation  of 
Pasquill  Category  A  conditions  in  MWT 
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Figure  9  Operating  range  for  CBL  simulations  of 
Pasquill  Category  B  conditions  in  MWT 
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Figure  10  Operating  range  for  CBL  simulation  of 
Pasquill  Category  C  conditions  in  MWT 


Figure  11  Combinations  of  surface  heat  flux  and 
wall  mass  transfer  required  to  produce 
specific  prototype  mixing  layer  growth 
rates 
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Figure  13  Dimensionless  velocity  profiles  versus 
surface  layer  position  for  finite  wall 
blowing  and  surface  heating  rates 


Figure  14  Dimensionless  temperature  profiles  versus 
surface  layer  position  for  finite  wall 
blowing  and  various  surface  heating  rates 
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Figure  17  wind-tunnel  subsidence  requirements  to 
simulate  the  Minnesota  CBL  experiments 
(Kaimal  et  al.,  1976)  for  various  model 
mixing  layer  velocities 


Figure  18  Wind-tunnel  wall  mass  transfer  rates 
required  to  simulate  the  Minnesota  CBL 
experiments  (Kaimal  et  al . ,  1976)  for 

various  model  mixing  layer  velocities 
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Figure  19  Calculated  velocity  profiles  anticipated 
in  the  surface  layer  when  simulating 
Minnesota  Case  2A1  {Kaimal  et  ai . ,  1976) 
by  augmented  wall  transport  methods 
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